Symmetry-breaking self-consistent quantum many-body structure of high-lying 
macroscopic self-trapped and superposition states of the gaseous double-well BEC 



o 
o 

(N 
-i— > 

o 
O 

m 

(N 



5h 



c 

o 
o 



> ■ 

o : 
\o . 
o ■ 

: 

o . 
-i— > ■ 

c3 ■ 



T3 

o 
o 



X 



David J. MasiellcO and William P. Reinhardf] 

Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, USA 

(Dated: February 6, 2008) 

The many-body structure of high-lying excited states, including macroscopic quantum superpo- 
sitions, of the gaseous double-well BEC is presented within the context of a multiconfigurational 
bosonic self-consistent field theory based upon underlying symmetry-broken one-body wave func- 
tions. To better understand our initial results, a model is constructed in the extreme Fock state 
limit, in which macroscopic quantum self-trapped and superposition states emerge in the many-body 
spectrum, striking a delicate balance between the degree of symmetry breaking, the effects of the 
condensate's mean field, and that of atomic correlation. It is found that, in general, the superpo- 
sition state lies energetically below its related self-trapped counterpart. Furthermore, noticeably 
different spatial density profiles are associated with each type of excited state. 
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Bose-Einstein condensates (BECs) of identical atoms 
confined to external double-well trapping potentials are 
actively being investigated Q, 0, EIEI • Remarkably sen- 
sitive matter- wave interferometry yj, based upon a co- 
herently split gaseous BEC source, represents just one 
potentially technologically useful experimental tool al- 
ready demonstrated in the double-well condensate. Re- 
cent experiments have also achieved exotic macroscopic 
quantum self-locked (self-trapped) states Q and, in do- 
ing so, open the way for direct experimental study of the 
fundamental nature of macroscopic quantum superposi- 
tions 0] directly observed in the gaseous BEC while only 
more indirectly in more conventional areas of condensed 
matter experiment 0, Q ■ None of the above experiments 
have been modeled in any detail at a first principles level. 

Theoretical descriptions of the double-well conden- 
sate's fully correlated many-body ground and low-lying 
excited states have been explored within the restricted 
Fock state basis \N U N 2 ) = (b\) Nl (S|) N * |vac) /y/Ni\N 2 \, 
where Nk is restricted to two macroscopically occupied 
states, k = 1,2. Here, the one-body wave functions (or- 
bitals) Xk approximating the bosonic field operators as 
*&(x) = xi(x)&i + X2(x)&2 have been chosen, with vary- 
ing levels of approximation, from solutions of the one- 
body Schrodinger equation [8j , to solutions of the bosonic 
Hartree-Fock equations 0, where the necessity of inclu- 
sion of orbital mean-field effects was established, as well 
as with solutions that have reached self-consistency in 
the occupation numbers n v k of the ^th quantum state 
[Icj . The roots of each of these approaches find their ori- 
gin in the quantum chemistry of many-electron atomic 
and molecular systems 1 1 . 

It is the purpose of this Letter to provide a novel 
theoretical first principles modeling of these high-lying 
excited states of the gaseous double-well BEC within 
the framework of an extended multiconfigurational 
bosonic self-consistent field (MCBSCF) theory that 
allows for symmetry breaking in both the many-body 



wave function and in the underlying orbitals \k, and 
includes a complete and consistent description of the 
effects of the condensate's mean field and atomic 
correlation in the restricted Fock space. In this the- 
ory, solutions of the many-body Schrodinger equation 
are constructed in which the orbitals underlying the 
TV-boson Hamiltonian H = J * t (x)/i(x)^(x)d 3 a; + 
(1/2) / *t( x )#t( x ')v( x ,x / )*(x')*(x)d 3 xd 3 a; / , have 
reached self-consistency with the many-body density. 
Within our restricted Fock space, the most general state 
of the system is a superposition of multiple configu- 
rations of identical bosonic atoms distributed, in all 
possible ways, between two condensate states at zero 
temperature. It is of the form 
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9«;{0'- i a}) = > C" Nl [C; a]\N lt N 2 ; {C; a}), (1) 



Ni=0 



where the nonlinear dependence of the orbitals (and the 
many-body state itself) upon the z4;h eigenvector of the 
Hamiltonian C v = {Cj^ : N x = 0, . . . , N} and the 
orbital symmetry-breaking parameter a are indicated 
in brackets. The novelty of this ansatz is that it ex- 
tends the work of Ref. to variationally include 
the effects of spatial symmetry-breaking, through a, in 
both the underlying orbital basis and in the macroscopic 
many-body wave function. The particular Fock basis 
state |JVi, N 2 ; {C v , a}) may represent, e.g., asymmetry- 
broken configuration with N\ atoms in the left well and 
N 2 atoms in the right well of the double-well trapping 
potential or perhaps a configuration with N± and N 2 
atoms in a symmetric and antisymmetric pair of states 
that enjoy the full symmetry of the Hamiltonian |l2l |. 
The specific nature of this basis and its connection to a 
will be clarified in the following. For a particular excited 
state v = 0, . . . , N, self-consistency is measured through 
the occupation numbers n v k , which are the eigenvalues 
of the one-body reduced density matrix 7 I/ (x,x / ), i.e., 
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J 7 !y (x, x')xfe(x')<i 3 a;' = n k Xk( x )- Both the number of of the fcth Fock state sum to the total number of atoms 
condensed atoms N k in and the occupation numbers n k N. By solving the coupled MCBSCF orbital equations 
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+ [/»(x)t& + v fcfc (x)i^ fcM + v«(x)rj; HI + v«(x)i^ HI + v ifc (x)r^„ fe ]^(x) 



for fe =^ i = 1,2, and then using these orbitals 
to construct the many-body Hamiltonian H, subsequent 
full diagonalization of H in Q yields many-body solu- 
tions 1^^; {C v , a}} that may, by iteration, reach self- 
consistency in occupation number with the orbitals un- 
derlying the vth quantum state. In Eq. @, h(x) = 
(— ft 2 /2m)V 2 + V cx t(x) is the one-body Hamiltonian, 
Vfcz(x) = / xSb(x')V"(x, x')xz(x')d 3 a;' are one-body ma- 
trix elements of the atom-atom interaction potential 
V(x,xM = (47rfi 2 a/m)(5(x— x') in the contact approxima- 
tion [13J with repulsive s-wave scattering length a, and 
j kl and T klmn (k,l,m,n — 1,2) are Fock space matrix 
elements of the one- and two-body reduced density ma- 
trices 7 w (x,x'J = (^(x^fx'^and P'(x,x';y,y') = 
(#t(x)*t( x ')#(y)$(y')) I/ 0,13 taken with respect to 
the MCBSCF state QJ. These matrix elements are func- 
tions of the expansion coefficients C u appearing in Eq. 

and, thus, are not specified but rather are connected 
to H through its eigenvectors. The field operators ^(x) 
and ^(x') an( i basic creation and annihilation operators 
bk and b\ satisfy the usual bosonic commutation rela- 
tions. 

Before discussing the MCBSCF theory of macro- 
scopic quantum self-trapped and superposition states 
of the BEC, it is instructive to consider a simplified 
model within multiconfigurational bosonic Hartree-Fock 
(MCBHF) theory @. By omitting the off-diagonal cou- 
plings in only the orbital sector of the MCBSCF theory, 
solutions of the resulting bosonic Hartree-Fock (BHF) 
equations are used to construct the many-body Hamilto- 
nian. Subsequent full diagonalization of H in Q yields 
solutions that may, similarly, achieve self-consistency 
through the diagonal parts of the density, or through 
application of the Hylleraas-Undheim theorem 0|. It 
is the latter approach that we employ in the following 
model while we strive for self-consistency in the general 
theory. 

Approximate implementation of restricted Fock space 
model. The diagonal terms in Eq. (0) form the general- 
ized BHF equations 

[ h lkk + VkkTkkkk + Viir%i kl \xk + VikFkiikXt = VkkXklkk 

(3) 

for N identical bosons where, again, k ^ I = 1,2. We 
point out that the matrix elements of 7^ and T u ap- 
pearing in Eq. (j3J) are not specified by BHF mean-field 
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theory, but rather find their origin in the full many-body 
theory where the state of the system may, in general, 
be a linear combination of several configurations. This 
fact will allow us, in the next paragraph, to derive sep- 
arate mean-field equations that are specifically tailored 
to model both types of aforementioned high-lying BEC 
states. 

In order to clarify and organize our discussion 
of macroscopic quantum self-trapped and superposi- 
tion states, we present a model that encapsulates 
the minimal physics necessary to describe their ba- 
sic properties. In the first case (a), all atoms 
are condensed into a single Fock state; the expan- 
sion coefficients for the z/th excited state are of 
the self-trapped form C v = [0, ... ,0, 1)^,0, ... ,0]. 
In the second case (b), the expansion coefficients 
mix two configurations with equal weight, i.e., 
C» = [0,...,0,(1/V2)^,0,...,0,±(1/V2)^ 2 ,0,...,0]. 
The resulting double-configurational many-body BEC 
state, which is the simplest multiconfigurational state 
imaginable, is a superposition involving a macroscopic 
number of atoms. We choose to consider these extremal 
model cases because MCBHF theory does not contain the 
full off-diagonal coupling of 7^ and Y w to the orbitals \k 
in the underlying (diagonal) BHF equations. However, 
in the extreme Fock state limit of (a) and (b), or even 
in slightly less extreme situations, the off-diagonal cou- 
plings contribute nothing or, at most, a small effect. All 
that remains is the diagonal coupling of the expansion 
coefficients to the orbitals, which we treat in this section 
exactly per manum by deriving specific BHF equations 
corresponding to each of the above situations. Such equa- 
tions will be explicated in the following. 

Corresponding to case (a), where the many-body state 
of the BEC is single-configurational, the diagonal Fock 
space matrix elements of the one- and two-body re- 
duced density matrices are j kk = N k , T kkkk = N k {N k - 
1), and T% lM = T" kllk = N k Ni, while in the double- 
configurational case (b) they are j kk = N/2, T kkkk = 
[N k (N k ~ l) +N l (N l - l)]/2, and = Tfa = 

NkNi. Substitution of these matrix elements into the 
general BHF equations © results in two sets of equa- 
tions that take into account the diagonal part of the 
coupling of the model states' expansion coefficients C v 
to the orbitals. The familiar BHF equations for iden- 



tical bosons |Tjj, llSl emerge in case (a), while in case 
(b), this coupling accounts for the diagonal interaction 
of the two configurations {l/\/2)\N\, iV 2 ; {C v ; a}) and 
±(1/^2)1^2, Ni,{C v ; a}} of the superposition state den- 
sity and results in a new pair of coupled mean-field equa- 
tions for superposition states. We solve the resulting sets 
of equations, mutatis mutandis, usin g a pseudospectral 
grid method in quasi-one-dimension 19] with standard 
relaxation methods (see Ref. @ for details) where, in 
addition, the orbitals are constrained to be of the form 

Xi (x) = sfaXi (x) + VI - «xf (x) 

(4) 

X2 (x) = VI -ax£(x) + V^f (x) 

with Xt an d Xk being the left and right half of Xk taken 
with respect to the trap center, and < a < 1 . This func- 
tional form allows for the orbital basis to be continuously 
fractionated from delocalized symmetric and antisym- 
metric pairs (a = 1/2) to symmetry-broken orbitals that 
are approximately localized in the left and right parts 
of the double- well trapping potential (a or a rj 1). 
One could imagine this variational parameter as being 
independent from the occupation number and perform- 
ing a second variation on n v k , however, for the purpose of 
our MCBHF model we further choose a to be consistent 
with n v k ; i.e., we take a = n±/N, which equals N±/N in 
this extreme Fock state limit. By varying a and then 
applying the Hylleraas-Undheim theorem |9j to each of 
the TV + 1 model MCBHF states \^ N ; {C";a = n y k /N}), 
we empirically find that this choice implies approximate 
self-consistency, by construction, for both low-lying and 
high-lying BEC states. This freedom of fractionation is 
motivated by the fact that the ground state and excited 
self-trapped and superposition state energies each find 
their minima at different values of a, or equivalently, in 
different orbital bases; the ground state and low-lying 
excited state energies are minimized in the delocalized 
(anti) symmetric basis, while the high-lying excited state 
energies reach their minima in the symmetry-broken lcft- 
and right-localized basis. It is well known that these 
bases are not unitarily equivalent in the mean-field ap- 
proximation for bosons |18j . 

Full diagonalization of H in the MCBHF states 
\^ N ;{C l/ ;a — n k /N}} corresponding to cases (a) and 
(b) results in the energy spectra displayed in Fig. Q] 
as a function of a = n\/N. In this, as in all other re- 
sults shown here, N = 100 and the scattering length a 
is chosen so that aN is representative of that value for 
typical laboratory BECs. The optimal MCBHF occu- 
pation numbers for the ground state and each excited 
self-trapped and superposition state, as variationally op- 
timized by the Hylleraas-Undheim theorem, are labeled 
with a bullet. Note that, in this approximation, the oc- 
cupation numbers appearing in the BHF equations may 
differ from those of the many-body calculation since off- 
diagonal couplings arise in the latter. For comparison 
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FIG. 1: (Color online) MCBHF energy E" as a function of 
fractionation parameter a — n"/N corresponding to macro- 
scopic quantum self-trapped states [case (a)] on the right, and 
macroscopic quantum superposition states [case (b)] on the 
left. In both cases, bullets lie at the optimal occupations 
for each excited state v as variationally optimized by the 
Hylleraas-Undheim theorem. For each self-trapped state v, 
the optimal fractionation a is determined by following curve- 
crossings. The total number of atoms TV = 100 and the details 
of trap geometry and appropriately scaled atomic interaction 
are similar to what is specified in Ref. @. Only the right 
half of _E"(0.5 < a < 1) for case (a) and the far left part of 
E"(0 < a < 0.15) for case (b) are displayed. 

with later MCBSCF results we display in Fig. [21 two 
pairs of BHF orbitals corresponding to the v = 90 self- 
trapped and superposition states. We point out that each 
of the high-lying superposition states have energies that 
lie lower than the corresponding self-trapped state of the 
same excitation, as can been seen in Fig. ^ For ex- 
ample, the v = 90 MCBHF self-trapped state, which is 
optimized at a = 0.88, has an energy of 261.373 Tilo while 
the equivalent superposition state, which is optimized at 
a = 0.055, has an energy of 256.938 hw. Due to the non- 
linear dependence of the energy upon the orbitals, this 
difference is a surprising new result. The spatial densities 

^(x) = 7 ,/ (x,x) = Xfc(x)7wXi(x) (5) 

H=l,2 

corresponding to these states as well as the MCBHF 
ground state, which is optimized at a = 0.5, are dis- 
played in the rightmost panel of Fig. |21 

Full implementation of restricted Fock space model. 
Having discussed both self-trapped and superposition 
states of the gaseous double- well BEC within the frame- 
work of our diagonal model, we now study these same ex- 
cited states within the more general theory that includes 
the off-diagonal couplings to the many-body densities 7" 
and r" in the orbital equations [see Eq. J2J], and, in 
addition, enjoys a variational degree of symmetry break- 
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FIG. 2: (Color online) MCBHF orbitals and densities corresponding to the ^ = 90 self-trapped state, optimized at a — 0.88, 
and the associated superposition state, optimized at a = 0.055. The self-trapped orbitals [case (a)] are displayed on the left 
and the superposition state orbitals [case (b)] are displayed in the center. In both cases, the nodeless orbital is xi and the 
nodal orbital is X2- They appear in the symmetric double well Vext at their appropriate chemical potentials. The far right panel 
presents the densities corresponding to each of these states as well as to the condensate ground state. 



ing. As before, the solutions of these MCBSCF equa- 
tions, which are parametrized by a, are obtained with 
standard relaxation methods. However, in this case, a 
is not linked to the occupation number; it is an inde- 
pendent variational parameter that controls the degree 
of spatial symmetry breaking in the underlying orbital 
basis [see Eq. (@J] and measures the fraction of each or- 
bital's L 2 -norm on the left- and right-hand sides of V ex t 
taken with respect to the trap center. Rather than ap- 
plying the Hylleraas-Undheim theorem to the resulting 
energy eigenvalues of H in order to determine the opti- 
mal occupation number for each MCBSCF self-trapped 
or superposition state \^ , {C v ;a}), we instead reach 
self-consistency in the eigenvalues n v k of the one-body re- 
duced density matrix 7"(x, x') separately for each v. The 
Hylleraas-Undheim theorem is then used to variationally 
optimize the MCBSCF energy E" as a function of frac- 
tionation a. 

Since the full coupling of and Y v to \k is included, 
there is no need to work within a model as we did pre- 
viously. The expansion coefficients C" may now involve 
many configurations in both Fock space and orbital sec- 
tors and will, by iterating until self-consistency, displace 
and distort the underlying MCBSCF orbitals through 
the many-body densities and repulsive coupling constants 
appearing in Eq. (0). In fact, self-trapped states owe 
their stability, in part, to an interference between neigh- 
boring configurations. The orbitals corresponding to 
the v = 90 macroscopic quantum self-trapped state, 
optimized at n\ — 11, a = 0.06, and the associated 
macroscopic quantum superposition state, optimized at 
n v k = 50, a = 0.035, are displayed in Fig. |3J The right- 
most panel presents the densities p v associated with with 
each of these states and, in addition, to that of the v = 
condensate ground state. In contradistinction to the or- 
bitals and densities displayed in Fig. it is clear that 
the MCBSCF based orbitals fill up the entire trap at a 
given chemical potential. Furthermore, each of the asso- 



ciated densities have noticeably different spatial profiles. 

In conclusion, we have developed a novel theoretical 
method, that is easily extendable to modeling actual ex- 
periments |l9j| . for characterizing the high- lying excited 
states of the gaseous double-well BEC that completely 
and consistently accounts for the inclusion of atomic cor- 
relation and condensate mean-field effects within a suf- 
ficiently general Fock space that is flexible enough to 
describe both macroscopic symmetry breaking and one- 
body orbital symmetry breaking. Many-body properties 
of both macroscopic quantum self-trapped and superpo- 
sition states of the BEC have been illustrated and differ- 
ences in their energies and densities have been pointed 
out. Our initial results suggest that Thomas-Fermi the- 
ory alone provides a good zeroth order approximation to 
the ground and high-lying excited state densities of the 
gaseous double- well BEC, suggesting the possibility of a 
hitherto unexpected simplification in the large N limit. 
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